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Integrability in Yang-Mills theory on the light cone beyond leading order 
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The one-loop dilatation operator in Yang-Mills theory possesses a hidden integrability symmetry 
in the sector of maximal helicity Wilson operators. We calculate two-loop corrections to the di- 
latation operator and demonstrate that while integrability is broken for matter in the fundamental 
representation of the SU{3) gauge group, for the adjoint SU{Nc) matter it survives the conformal 
symmetry breaking and persists in supersymmetric Af — 1, M = 2 and A/" = 4 Yang-Mills theories. 
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Four-dimensional Yang-Mills (YM) theory exhibits a 
new symmetry which manifests itself on the quantum 
level through integrability property of the dilatation op- 
erator in the sector of the maximal helicity Wilson op- 
erators in the multi-color limit Q. In supersymmetric 
Yang-Mills (SYM) theories, integrability is promoted to 
a larger class of operators and ultimately to all Wilson 
operators in the maximally supersymmetric TV = 4 the- 
ory 13] ■ It is expected that the A/" = 4 dilatation oper- 
ator has to be integrable to all orders in 't Hooft cou- 
pling since the TV = 4 SYM is dual to the superstring 
theory on AdSsxS^ background 0, whose world-sheet 
sigma-model possesses an infinite number of integrals of 
motion both on classical 0| and quantum levels^- 
far integrability in YM theory was unequivocally estab- 
lished at leading order in coupling and it was argued to 
hold in certain closed subsectors of the TV = 4 SYM the- 
ory in higher loops Q. In distinction with the TV = 4 
model, the conformal symmetry of YM theory is broken 
by quantum corrections and the conformal anomaly af- 
fects the dilatation operator starting from two loops. A 
natural question arises whether integrability in the sec- 
tor of the maximal helicity operators carries on to higher 
orders in YM theory and its supersymmetric extensions. 

In this Letter, we report on a calculation of the two- 
loop dilatation operator in the sector of three-quark 
(baryon) operators of the maximal helicity-3/2. Inte- 
grability is a genuine symmetry of YM theories and the 
choice of the sector is driven solely by the simplicity of 
the computations involved. For our purposes we con- 
sider separately the case of quarks in the fundamental 
S't/(3) and adjoint SU{Nc) representations. The former 
corresponds to QCD while the latter is relevant to its su- 
persymmetric extension. The corresponding three-quark 
light-cone operators are defined as 

Bf (zi , Z2 , ^3) = {zi)q] {Z2)ql{z3), 

Ba(zi, Z2, Z3)=tr [q^Zi)q^Z2)qHz3)] , (1) 

where q'^ (z) = ^(1 — 75)(7(znp) is the helicity-1/2 fermion 
"living" on the light-cone (n^ = 0) and, in the second 
relation, q — q°-t°- with the SU{Nc) generators in the 
fundamental representation. The quarks in |^ may have 
different flavors and their total number is Nf . It is tacitly 



assumed that the gauge invariance is restored in |^ by 
inserting appropriate Wilson lines between quark fields. 
Later on we shall adopt the light-like axial gauge n-A = 
in which the gauge links reduce to the unit matrix. The 
light-cone operators are generating functions of local 
Wilson operators. The latter can be obtained by Taylor 
expanding M{zi) in the light-cone separations z^. The 
operator M{(zi) has a direct phenomenological interest in 
QCD as its matrix element defines the A'^/'^— distribution 
amplitude .7]. For Nf = 1, the operator Mi,(zi) is a 
component of a supermultiplet in the TV = 1 SYM theory. 
Its remaining components are obtained from Ba(zi) by 
replacing q^ one by one with helicity-up gluon fields, thus 
covering the maximal-helicity three-gluon operator. 

The goal of this study is to elucidate the symmetry 
properties of the dilatation operator governing the scale 
dependence of the operators to two-loop order. A 
unique feature of the operators ^ is that they mix only 
with themselves and obey the Callan-Symanzik equation 



i{z,) = -[H-B](z,) 



(2) 



Here H is an integral operator acting on light-cone coordi- 
nates of quark fields and defining a representation of the 
dilatation operator D in the space spanned by light-cone 
operators, z[D,B(z,)] = (M + Mq)M{z,) with = 3/2 the 
quark canonical dimension. Its perturbative expansion 

H = AH<'"+A2h"' + C»(A3), (3) 

goes in A, depending on the quark representation in 



4 g'^ 
^^^8^3' ^^^8^^^=' 



(4) 



and gauge coupling g. Notice that the leading order evo- 
lution kernel H'"' is universal (see below) while the 
two-loop kernel H*^' is different for Bf (z;) and Ba(zi). 

To solve (0) one has to examine the Schrodinger equa- 
tion for the evolution kernel ||2Jl 

[H • *5](zi, Z2, Zs) = 79(A) '^qizi,Z2, Z3) , (5) 

with eigenvalues 7g(A) having a perturbative expansion 
in the coupling constants Q similar to Eq. and 

*,(z,) = vl/<")(z,) + \^'^\z,) + 0{X') (6) 
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FIG. 1: Two-loop corrections to the two-particle (o — /) and three-particle (g, /i) evolution kernels in II12|I . 



being homogeneous polynomials of degree = 0, 1, . . .. 
Having solved one can expand nonlocal operators 
B(zi) over local Wilson operators 0^(0) built from three 
quark fields and A'' covariant derivatives acting on them 



*5(Z1,Z2,Z3)0),(0). 



(7) 



To one-loop order, these operators have autonomous 
scale dependence but this property is lost in two loops 
due to dependence of ^^(zi) on the coupling constant, 
Eq. Integrability arises as a hidden symmetry of the 
Schrodinger equation jSjl through the existence of a new 
quantum number q. 

To one-loop order, the kernel has a pair-wise form 
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with the two-particle kernel Elj^.'' displacing j"^ and k 
quarks along the light-cone in the direction of each other 



(8) 



,th 



[H<"2>.B](zi,z2,;23)= / —a 
/o a 



1 da 



2B(zi,Z2,23) (9) 



l{azi + az2,Z2, Z3) ~ B(zi, azi -f az2, Z3) 



where a = 1 — a and jq = 1 is the conformal spin of the 
quark. The one-loop dilatation operator ||SJ| possesses 
two conserved charges 4"', g^"' such that [H<°', ] = 0, 

=i?2+ ^23+^31, ^^[Ll2,Ll,], (10) 

where L^f. is a two-particle Casimir operator of the 
"coUinear" SL{2) subgroup of the conformal group 



order, we calculate the two-loop dilatation operators © 
in the light-like axial gauge n ■ A = 0. In this gauge, the 
relevant Feynman diagrams are shown in Fig. ^ Their 
contribution to the two evolution kernels H^"^' and H*^^' 
only differs by the accompanying color factors. As a con- 
sequence, the two kernels are different but have the same 
general form 



= Jnii2 



1123 



1131 



II123 



(12) 



Here Mjk is the two-particle kernel computed from the 
graphs in Fig.n(a) - (/) and H123 is the three-particle ir- 
reducible kernel defined by the diagram {g) and its cyclic 
permutations. The diagram (/i) vanishes for both op- 
erators in 1^. Notice that the nonplanar diagram (b) 
contributes to Bf (z^) and vanishes for Ba(zi). The quark- 
leg renormalization (not displayed) contribute to the c- 
number constant F'^'. 

By analogy with one- loop, Eq. (jS)), we define the two- 
loop Hamiltonians Mjk and M123 in such a way that they 
annihilate the local baryon operator Ba,f (zi — 0) so that 
its anomalous dimension equals [H — 7o(A)]B(zi = 0) = 
with 70(A) = A f F<i' + 0(A3). To save space we do 
not present the explicit expression for F'"^'. It defines the 
ground state energy for the Hamiltonian H, Eq. |(SJ, and 
is not relevant for the discussion of integrability. 

The diagram in Fig. ^ (/) involves a quark loop. Its 
contribution is proportional to the number of flavors Nf 
and is gauge invariant. This allows one to split the two- 
particle kernel Mj^ into the sum of gauge-invariant terms 



(A/ck): 



JT(l) 



(13) 
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(11) 



with Zjk = Zj — Zk- Complete integrability implies that 
the one-loop dilatation operator is a function of the con- 
served charges H*"' = H<'"(4°', g'"'). So that its eigen- 
functions diagonalize both of them simultaneously, e.g., 
g^'"«'^'"(-z») = Q'^q\z,). Remarkably, H"" can be mapped 
into the SL{2,R) Heisenberg magnet of spin— jg and 
length equal to the number of quarks 0. As a result, the 
Schrodinger equation Q is completely integrable to one- 
loop order and is solved exactly using the Bethe Ansatz. 

The one-loop integrability holds irrespective of the 
SU{Nc) representation of the quark fields in To ver- 
ify whether this symmetry is preserved beyond leading 



where the second one is proportional to the one-loop 
function f3o = lliVc/3 - 4rniV//3 in the underlying YM 
theory with Tf = 1/2, Cf = 2/3 and Ta = Ca Nc/2. 

Starting from two loops the anomalous dimensions re- 
ceive conformal symmetry breaking corrections due to 
nonzero /^-function H13I) and, in addition, they depend 
on the scheme used to subtract ultraviolet divergences. 
In our calculations we applied the dimensional regular- 
ization with d = A — 2e and subtracted poles in 1/e in 
the MS— scheme. In general, supersymmetry is broken 
in the dimensional regularization and is preserved in the 
dimensional reduction. The diagrams in Fig.^give iden- 
tical results for both procedures so that supersymmetry is 
preserved. The difference arises only for two-loop quark 



self-energy, which contributes to F'"^'. The resulting ex- 
pressions for the two-loop Hamiltonians Wjk and IHI123, 
Eq. 1)12(1 . are cumbersome and will be given elsewhere. 

The dilatation operator H, Eq. J^Jl, is invariant under 
the cyclic permutation of quarks P and the interchange of 
any pair of quarks Fjk- Since [P, Pjk] 7^ 0, the eigenfunc- 
tions ^!q(zi) cannot diagonalize both operators simulta- 
neously. Choosing P, the eigenstates can be classified 
with respect to the quasimomentum 

P«'q(zi,Z2,Z3) = *g(2:2,Z3,Zl) = e'''-'^g{zi,Z2,Z3) 

(14) 

with 6q = 27rn/3 and n — —1,0,1. The parity symme- 
try of the Hamiltonian, [Pi2,H] = 0, combined with the 
identity PP12 = Pi2P^ implies that, independently of 
the existence of the charge 93, the eigenvalues of H with 
non-zero quasimomentum 9q are double degenerated, 



[H-7,(A)]vI/±(z,)=0, vff = i(l±Pi2)*,(z.), (15) 

whereas for 9^ = the eigenvalues are not necessarily 
degenerated. The existence of the conserved charge (73 — 
93°' -|- A ^3^' -|- . . . symmetric with respect to P but antisym- 
metric with respect to P12, i.e., [P, (73] = {Pi2,93} = 0, 
extends the degeneracy property to the eigenstates with 
9g = and qs + 0, since Pi2H((73)Fi2 ^ H(-g3). 

At one-loop, integrability of the Schrodinger equation 
(O implies double degeneracy of all eigenvalues of IHI<°' , 
except those with (73 = 0. To determine two-loop cor- 
rections to the eigenspectrum of H, Eqs. H12|) . we apply 
the conventional 'degenerate perturbation theory'. Using 
the basis of one-loop definite-parity eigenstates ^'^'''^(zi), 
Eq. (|15|l . normalized with respect to the 5^(2; R) invari- 



ant scalar product as ) 

= |H<''|*<'"±}, 



bqqi , one finds 
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(16) 
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Going over through lengthy calculations we observed 
that, in agreement with our expectations, the eigenstates 
with nonzero quasimomentum 9q ^ Q are double degener- 
ated to two loops for both operators in (QJ . For the eigen- 
states with 9q = Q the situation is quite different. We 
found that the desired pairing of eigenvalues occurs only 
for the operator Ba(zi) while it is lifted for the Bf(zi). 

We remind that the eigenstates ^'^(zi) are homoge- 
neous polynomials of degree N = 0, 1, ... and the total 
number of eigenvalues equals + 1. The first non-trial 
example arises for iV = 3 and demonstrates the main 
trend, recurring for higher N . At TV = 3, for quarks 
in the adjoint representation one finds two pairs of the 
eigenvalues A7=^ = 7^ (A) - 70(A), 



A7f(A) = Aa4 + A2 



29 

24 



(17) 



with A7j^(A) and A7j^(A) corresponding to the eigen- 
states with 9q — Q and 9q = ±27r/3, respectively, and 
/3o/ca — §(11 ~ 2A^/). For quarks in the fundamental 
representation, the two eigenstates with 9q = Q have dif- 
ferent eigenvalues starting from two loops, A7j''' 7^ A7j~, 
and we get instead 



A7j+(A) = Af4 + Af' 
A7f (A) = Af4 + A^ 
A75(A) = Af^ + A2 
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(18) 
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Nf. For Nf — 1 the exact spectrum 



with /3o/cf 

of eigenvalues A7(A) = A7]^' + A^7]y' for the conformal 
spin < A^ < 10 is displayed in Fig. [21 It clearly demon- 
strates the lifting of the degeneracy and proliferation of 
eigenvalues for the baryon operator Bf(zi) compared to 
Ba(zi) where the pairing persists for all eigenstates. 

We recall that the conformal symmetry is broken for 
/3o 7^ 0. One immediately sees from (I17|l and H18|l that 
though the degeneracy is lifted for the operator Bf even 
in the conformal limit /3o = 0, the /3o-terms themselves 
do preserve the "pairing" of eigenvalues. We have to em- 
phasize that for quarks in the adjoint representation, the 
exact expressions for the anomalous dimensions (|17ll co- 
incide with those in the multi-color limit. At the same 
time, for quarks in the fundamental, the baryon opera- 
tor B{{zi) exists only for Ac = 3, so that the large- A'c 
counting is not applicable. The non-planar diagramn(b) 
contributes to 1(18(1 and partially leads to breaking of in- 
tegrability in two loops. 

Remarkably enough, the degeneracy of the eigenstates 
of the two-loop dilatation operator H survives the con- 
formal symmetry breaking. To understand the phe- 
nomenon, one considers the dilatation operator H in 
the large /3o— limit, XPo = fixed and Pq — > 00. In 
this limit, the {n + 1)— loop Hamiltonian H is domi- 
nated by the contribution of A(A/3o)"— terms, which can 
be resummed. In two loops, A^/3o— correction to l(T^ 
comes from the second term in (II 3() involving the kernel 
Being combined with the one-loop kernel, the sum 

operator as 



jk ■ 



is determined by the same integral 
!, Eq. ©, 



provided that one substitutes 

, /3o.g' /5 



167r2 



ln( 



(19) 



Adding higher-loop A(A/3o)"— corrections, one finds that 
the In a— terms in the right-hand side of ((19(1 exponenti- 
ate and additively renormalize the bare quark conformal 
spin as jq — > jq-t-/3o5^/(327r^) 's']. As a result, the all-loop 
dilatation operator H, Eqs. (O, is given in the large— /3o 
limit by the one-loop expression (|SJ) with the quark con- 
formal spin jq = 1 replaced by its 'renormalized' value 
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FIG. 2: Spectrum of anomalous dimensions (see text) at one loop (a), and two loops for quarks in the adjoint (b) and 
fundamental (c) representations. In the latter case, more states for the same conformal spin N arise due to lifting of the 
degeneracy. 



with ip{x) = Qr(i-''x)i^H^+x) ' = '^f and Aa for the op- 
erators M[{zi) and Me,{zi), respectively. The operator 
]g[(/3o^oo) inherits integrability of the one-loop dilatation 
operator IHI<°' and it can be mapped into the Hamiltonian 
of the Heisenberg magnet of spin jg = 1 + /3og^/(327r^). 
It possesses two conserved charges q2°^°°^ and q^f°'^°°^ 
which are given by the same expressions as before, 
Eq. ifTUI) . where L^^. is given by ((TT|) with = 1 + 
/3o5^/(327r^). Thus, the /Sg— terms preserve integrabil- 
ity of the dilatation operator for both operators in 
and do not lift the degeneracy of its spectrum. 

For Po = 0, integrability holds to two loops only for 
adjoint quarks. Indeed, introducing the operator Z = 
(ll-FPi2)Z+/2-h(l-Pi2)Z-/2 with Z± defined in (Cnj, 
it is straightforward to verify that the charges 

gn=C+Aa[Z,C]+0(AD, (" = 2,3) (21) 

satisfy [92,93] = [92, H] = [93, H] = + 0{Xl) provided 
that the spectrum of H is degenerate, 7^^*"'' = 7^^'". 
Notice that the two-loop dilatation operator H in the 
MS— scheme does not respect the conformal symmetry 
even for f3{g) = 0, [H, 93°'] 0- The reason for this 
is that the dimensional regularization inevitably breaks 
conformal symmetry since P{g) — {d — 4:)g/2 is non-zero 
in dimensions and it propagates into the anomaly as 
c? ^ 4. The anomaly can be cured by a finite renormal- 
ization of the B— operator, 8(7(2^) = [U-B](zi), leading 
to the dilatation operator, Wu = U-H-U-^ [3. This 
does not affect the spectrum of the anomalous dimen- 
sions but modifies the conserved charges 92 , Q3 and the 
eigenfunctions 'i'q{zi). The operator U defines a scheme 
in which the conformal symmetry of the dilatation oper- 
ator is restored for (3{g) = 0. Such scheme is not unique 
and can be defined for our best convenience. Choosing 
U = 1 - AZ + 0(A2) with the same Z as in one 
can define a scheme in which the charges do not receive 
radiative corrections {qn)u ~ VqnV^^ — -I- ©(A^) 
while the two-loop correction to the dilatation operator 
Mu = Il<") + AH<)' -I- ©(A^) with H<}' = H<'' - [Z,H<'"] 
commutes with the one-loop kernel, []HI'°',]HIj)'] = 0. 



So far we discussed the three-quark operators As 
was already mentioned, at Nf = 1 the operator Ma{zi) is 
a component of the supermultiplet in the M — 1 SYM. 
Supersymmetry allows one to extend the two-loop inte- 
grability to remaining components of the supermultiplet 
including three-gluon operator of the maximal helicity 
Going over to TV = 2 and A/" = 4 SYM theories, one 
finds that the only modification one has to make is to add 
to the dilatation operator (|12(l the contribution of two- 
loop diagrams with, respectively, two and six real scalars 
inside the loops of two-particle kernels but not in the 
external lines since the maximal-helicity operators carry- 
ing the maximal i?— charge mix only among themselves. 
The scalars modify the constant r<^> in (|12|l and induce 
an additional contribution to the two-particle kernels Wjk 
proportional to H^°j! and H^^j! . Since these terms preserve 
integrability, the two-loop integrability found above is 
promoted to the TV = 2 and TV = 4 SYM theories. 
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